Piezoelectric control of the magnetic anisotropy via interface strain coupling in a 

composite multiferroic structure 
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We investigate theoretically the magnetic dynamics in a ferroelectric/ferromagnetic heterostruc- 
ture coupled via strain-mediated magnetoelectric interaction. We predict an electric field-induced 
magnetic switching in the plane perpendicular to the magneto-crystalline easy axis, and trace this 
effect back to the piezoelectric control of the magnetoelastic coupling. We also investigate the 
magnetic remanence and the electric coercivity. 
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PACS numbers: 75.85.+t,77.65.-j,85.70.Ec,75.70.Cn 

Controlling magnetism with an electric field is funda- 
mentally important and bears the potential for a wide 
range of applications as for instance in sensorics and mag- 
netoelectrically controlled spintronics devices with ultra 
low heat dissipation Multiferroics, i.e., materials 

that exhibit ferroic (magnetic, electric, and elastic) or- 
ders allow to tune the coupled ferroelectric (FE) 
and the ferromagnetic (FM) order parameters by exter- 
nal magnetic and electric fields, respectively. A major 
obstacle however is the relatively low strength of the mag- 
netoelectric coupling in bulk matter. A promising route 
to circumvent this problem are the appropriately synthe- 
sized composite FE/FM nano and multilayer structures 
that may serve as elements in quantitatively new multi- 
ferroic devices at room temperature [6-8]. One example 
that has been studied recently theoretically and experi- 
mentally is BaTi0 3 (BTO)/Fe [1-QJ]. Steering the mag- 
netization of Fe layer in BTO /Fe by an electric field was 
demonstrated experimentally in Refs. [1, [HI, [3 . Several 
mechanisms may underlay the magnetoelectric effect in 
composite multiferroic junctions, for example ch arg e re- 
arrangements [1, 0, Jl m l 3 , strain effects [13, HH, [3 Ell > 
and exchange-bias [20l - [22j have been studied as the key 
ingredients for the magnetoelectric coupling. In a pre- 
vious study we investigated theoretically the dynamic 
response to an applied external field in a multiferroic 
chain with a linear magnetoelectric coupling that results 
from the electrostatic screening at the FM/FE interface 
[l2l [23|]. An important finding is that for the material 
parameters corresponding to BTO/Fe the total electric 
polarization and the net magnetization are controllable 
by external magnetic and electric fields, respectively [24j . 

An outstanding important problem however is the 
question of how the strain, due to the interface lattice 
distortion in a FE/FM heterostructure, affects the mul- 
tiferroic dynamics. Experimentally the thin FM Fe film 
on BTO single-crystal was successfully realized, and the 
magnetic properties of Fe film are found to be strongly 
modified by the successive structural transitions of BTO 
substrate via interfacial magnetoelastic coupling in three 



different phases, viz., tetragonal (300 K), orthorhom- 
bic (230 k), and rhombohedral (150 K) BTO phases 
fiol [l3l IT5I [Hj]. All these experimental finding are 
very promising steps towards the manipulating magnetic 
anisotropy using interfacial strain. However, the ferro- 
electric control of magnetism by the structure phase tran- 
sition via the thermal cycle lacks high-speed read/write 
multiferroic dynamics as to their intrinsic relative long re- 
laxation time. In contrast, the thin BTO film possesses 
strong piezoelectricity in experiments [HI, H3- Epitaxial 
strain are found to enhance the spontaneous polariza- 
tions of the FE thin film. Furthermore, as shown by 
Ederer and Spaldin in Ref.j26j, the strain dependence 
of the polarization in thin FE film can be understood 
in terms of the piezoelectric and the elastic constants of 
the unstrained materials. In this letter, we present the- 
oretical results for the interface strain-based magnetic 
response to an oscillating external electric field based on 
the piezoelectricity, instead of the thermal phase tran- 
sition. The system under consideration is illustrated in 
FigfT] An cpitaxially grown multiferroic heterostructure 
that we treat in a coarse-grained manner as consisting of 
Nfm ferromagnetic (e.g. Fe) and A/fe ferroelectric (e.g. 
BaTiOa) sub-cells. Each cell is a cube of equal volume 
a 3 — 5 x 5 x 5 nm 3 , as in Ref . [24j . 

For the epitaxially strained FE films, the strain com- 
ponents satisfy thus 7/4 = 715 = uq = 0, u\ = ^ 0, and 
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FIG. 1. Schematic view of the ferroelectric/ferromagnetic het- 
erostructure. The BTO film is assumed [001] oriented with 
the polarization perpendicular to the interface (pointing ei- 
ther e z or —e z ). The magnetic moments are taken with all 
three cartesian components and the initial directions are cho- 
sen randomly. 
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the strain perpendicular to the interface is determined by 
the Poisson ratio n = —ui/u 3 , where i = 1, 6 denotes 
the strain tensor in Voigt notation. Up to a linear order, 
the change in the polarization P is given by the improper 
piezoelectric tensor, c a i = [23]. On the other hand, 
the magnetoelastic coupling at the interface is associated 
with an additional uniaxial anisotropy energy for the FM 
layer [Hj], 
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where A is the average magnetostriction coefficient, and 
<\> is the angle between the magnetization M and the 
direction of the stress er. With these consideration of the 
stress-strain effect, the total free energy density of the 
multifcrroic chain is 

Fmf = Ffe + Ffm- (2) 

The energy density of the FE subsystem is given by 

Ffe = Fq + Fddi (3) 

where the elastic Gibbs function Fq reads 

Fa = - a £ P 2 + p £ Pj + k £(P, - P,-!) 2 

i i i 
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Here c e g = 2c3i — 033/71 (note that the symmetry C31 = 
C32 has been exploited to derived the effective piezoelec- 
tric constant [26[ ) , C u is the stiffness coefficient of the FE 
film, and E(<) = (0, 0, Eg sin u>t) is an applied harmonic 
electric field. The long range FE dipolc-dipole interac- 
tion Fddi has the usual form 
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where epE is the FE permittivity, eo is the dielectric con- 
stant of vacuum, is the distance between P, and P^, 
and eij is the unit vector joining the two dipoles. 

For the ferromagnetic energy density the relation ap- 
plies 



Ffm — Fxc + Fmst + Fj 
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Fxc consists of the nearest-neighbor exchange interac- 
tion and the uniaxial magneto-crystalline anisotropy con- 
tributions 
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with the saturation magnetization M s . The magnetic 
dipolc-dipole interaction is 
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Fmmi = 7— / 
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where /xo is the magnetic permeability constant. 

For the FE subsystem, the material parameters are 
chosen as a = 2.77 x 10 7 Vm/C [pL /3 = 1.70 x 10 8 
Vm 5 /Cl_S, k = 1.0 x 10 8 Vm/C 13, and P s = 0.499 
C/m 2 [2J]. The improper piezoelectric constants are set 
as those of BTO c 3 i = 0.3 C/m 2 , c 33 = 6.7 C/m 2 , 
and the Poisson ratio n = 0.64. The average stiffness co- 
efficient C u of thin FE film is unknown and expected to 
quite smaller than that in the bulk [l8j . In the numerical 
simulation we used C u = 3.90 x 10 9 N/m 2 that is adjusted 
to reproduce the lattice mismatch (~ 1.4%) between the 
bcc Fe[001] film and the tetragonal BTO at room temper- 
ature [15j . Further material parameters concerning the 
FM layer are iron, i.e., A = 2.07 x 10~ 5 along Fe [100] 
2j|, A = 2.1 x 10" 11 J/m [Iljjfi = 4.8 x 10 4 J/m 3 
3(|, and M s = 1.71 x 10 6 A/m (3^. The induced stress 
acting on the FM body is given by er = —C\\{u\, U2,Us) 
with C\\ = 1.78 x 10 11 N/m 2 being the elastic constant 
of BTO at the interface [3l|. 

The multiferroic dynamics is studied by Monte Carlo 
simulations with the standard Metropolis algorithm [32[ 
and open boundary condition at T — 312 K for tetrag- 
onal BTO phase. The FE dipoles P,; are chosen to be 
uniformly distributed over (0, 0, — P s )-(0, 0, P s ) by a ran- 
dom function [33l [34[ . The three-dimensional magnetic 
moments M.; are updated coherently, i.e., at each trial 
Monte Carlo step a direction of new is limited within 
a cone around the initial spin direction 35[. The maxi- 
mum angle 9 ma x of the cone is determined by means of a 
feedback algorithm so that the number of accepted spin 
modifications is just half the total number of equilibrium 
configurations at temperature T and at time t = in 
the absence of the external electric field E(t) [36|. Af- 
ter the multiferroic equilibrium at T = 312 K is estab- 
lished, 9 max is fixed and a sinusoidal oscillating electric 
field along the chain, E z (i) = Eq sinqjt, is turned on and 
sweep through the multiferrroic chain [37h39| . One sweep 
constitutes one Monte Carlo step (MCS) , which is taken 
as the time unit. To reduce random errors, the data are 
collected and averaged for 400 full electric field cycles. 

Figf5]shows the hysteretic piezoresponse of the FE sub- 
system. The FE layer is poled and expands or shrinks in 
response to the electric field. Experimentally a similar 
piezoresponse has indeed been observed in BTO/LSMO 
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FIG. 2. FE polarization switching (left) and piezoresponse 
loop (right) driven by the harmonics electric field with period 
2-k/u = 200 MCS. Hysteresis curves are calculated for Nfe = 
20 and N FM = 4. 
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FIG. 3. Multiferroic response of the magnetic anisotropy to 
applied electric field. M a is the projection component of mag- 
netization in the direction of interface stress. The parameters 
are the same as in Fig(2] 



sample [14[. In BTO film, the dynamics of strain is di- 
rectly related to phonon modes of the crystal lattice. 
Their launching time in each cell is around ~ 1 ps, 
estimated with the sound velocity veto ~ 5000 m/s 
at the room temerpature [40j . which is much shorter 
compared to the oscillation period of the applied elec- 
tric field (HI, We have thus a linear piezoelec- 
tric dynamic coupling of ferroelectric and strain modes 
and time-independent Poisson ratio. This piezoresponsc 
in turn updates the induced interfacial stress, cr(t) = 
— CiiUi(t)(l, 1, — 1/n), resulting in a time-oscillating uni- 
axial magnetic anisotropy, consequently an anisotropic 
multiferroic behavior in the FM film (c.f. FigJ3]). The 
maximum of the stress anisotropy energy is estimated as 
Fmst = 1-7 x 10 5 J/m 3 , which gives rise to a pronounced 
magnetic response to the electric field in the normal- 
plane and along the direction of the stress <r even though 
the modification of the total magnetization is quite small. 
As Mi > 0, we have U3 < that implies a tensile strain to 
Fe along the chain. The stress axis becomes magnetically 
harder, which favors a small alignment of the magneti- 
zation along the stress. Whereas, u\ < indicates a 
compressive strain (Act > 0), that gives rise to an easy 
cr-axis and a large projection M a (c.f. Fig|3j Such an 
anisotropic magnetization behavior is q ualitatively con- 
sistent with the experiment [T3, EH, 1 1 Of ] . The change of 



in-plane magnetization reads AMy/Af, ss 11%, smaller 
than that (~ 33%) induced by the structure phase tran- 
sition |10j. However, an advantage of the piezoelectric 
control is that the magnetic anisotropy can be changed 
by a high-frequency oscillating electric field instead of the 
slow thermal cycle. 

It should be noted that the influence of the stress-strain 
effect on the magnetic anisotropy is typically more con- 
siderable than the surface magneto-crystalline anisotropy 
stemming from the charge redistribution at the FM/FE 
interface [1, EH EH EH Ell- And the anisotropy intro- 
duced by the interfacial spin-polarized screening charges 
is along the chain [IH [2J|, which results in a periodic 
variation of M x and M y but does not allow to form 
the hysteresis loop (24|. To further characterize the 
electrical switchability of the in-plane magnetization of 
the FM film, we study the frequency dependence and 
the finite-size effect on the magnetoelectric response, as 
shown Fig|4] The electric coercive fields for M x hystere- 
sis loop are plotted as a function of the frequency of the 
applied electric field. Clearly, the increase in the fre- 
quency results in a decrease in the coercivity field E c . 
For thicker FM layer, the interface strain effect becomes 
weaker and is eventually overwhelmed by the magneto- 
crystalline anisotropy, as marked by the fast decrease of 
the remanent magnetization M x with Nfm- 

In conclusion, we demonstrated theoretically that the 
interface strain effect can strongly influence the multifer- 
roic dynamics in FE/FM heterostructures. Through the 
piezoelectricity and the magnetostriction, the induced, 
electrically tunable magnetic anisotropy results in a well- 
developed in-plane magnetic hysteresis loop as a function 
of the electric field. This finding has a promising po- 
tential for the electric control of the complete magnetic 
switching whilst the total magnetization remains stable 
in presence of applied voltage. In order to get integrally 
insight into the multiferroic dynamics, the magnetoelec- 
tric coupling induced by the screening charges should 
be taken into account as well. For the strength of the 
magnetoelectric coupling in the range above 1 s/F [24j . 
two types of multiferroic interactions are comparable. 
The stress anisotropy energy prefers the in-plane mag- 
netization, whereas the surface spin-dependent screening 
charges favors to antiparallelly align the FE polarization 
and FM magnetization along the chain. Studies of the 
competition between them promises rich potential for fu- 
ture multiferroic devices and are currently under way. 
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FIG. 4. The electric coercivity (left) and magnetic remanence 
(right) for M x hysteresis loop of the FM/FE heterostructure 
with 20 FE cells. 
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